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We study the quasiparticle resonance states near strong impurities in the mixed state of a d-wave 
superconductor. These states give rise to zero-bias peaks in the local density of states, observed 
in scanning tunneling microscopy experiments. The field dependence of the peaks is obtained by 
averaging with respect to the spatially non-uniform Doppler shifts in the energy of excitations. The 
hybridization of the magnetic field-induced nodal quasiparticles with the impurity resonances results 
in the suppression and broadening of the zero-bias peaks. The height of the peaks is found to scale 
as {H c2 /H) 1/2 . The magnetic field response of the peaks is shown to be strongly dependent on the 
field orientation. 
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One of the striking features of high-T c superconductors 
(HTSC's) is that even a single impurity has a notable 
effect on the superconducting state, creating a sharp res- 
onance state in its vicinity, which becomes a true bound 
state of zero energy in the limit of strong impurity po- 
tential (the unitary limit) [3, E H & H, @ . These bound 
states manifest themselves as sharp peaks near zero bias 
in the energy dependence of the single-particle density of 
states (DoS), which have been observed in scanning tun- 
neling microscopy (STM) experiments on BSCCO com- 
pound HSUEi 

The bound states near non-magnetic impurities are 
specific to unconventional superconductors with an 
anisotropic order parameter and the gap nodes at the 
Fermi surface 11]. The zero width of these states in 
the unitary limit follows from the absence of hybridiza- 
tion with the extended quasiparticle states: HTSC com- 
pounds are d-wave superconductors, and the DoS, N (uj), 
of the extended states vanishes at small energy linearly 
in lu. However, this is no longer true if the superconduc- 
tor is placed in an external magnetic field H . The field- 
induced supercurrents act as pair breakers, filling the gap 
nodes and creating a finite density of bulk quasiparticles 
at low energies. These quasiparticles are responsible for 
a number of peculiar properties of the HTSC's, such as 
a non-analytical in H behavior of the electronic specific 
heat and the thermal conductivity in the mixed state at 
low temperatures (the Volovik effect) |l3 ]. One can ex- 
pect that the hybridization with the field-induced bulk 
quasiparticles gives rise to a finite width of the impu- 
rity bound states. A simpler version of this problem, 
applicable, for instance, to the case of H || ah in the 
Meissner state, was studied in Ref. 0, where a uniform 
supercurrent q s was considered and it was shown that 
the bound states indeed acquire a non-zero width which 
depends non-analytically on q s . In this paper, we study 
the effect of the magnetic field on the resonance states 
near unitary impurities in the mixed state for H || c. In 
the presence of vortices, the supercurrent is non-uniform, 
which qualitatively changes the magnetic field response 



of the impurity resonances. The unitary limit of scatter- 
ing is of particular interest because the most profound 
effects related to the impurity resonances have been ob- 
served in the vicinity of Zn impurities in BSCCO, which 
have the s-wave phase shift Sq close to ir/2 Q. 

Suppose we have a repulsive point-like impurity which 
is described by the potential U(r) = u5(r) (u > 0) in 
a two-dimensional d-wave superconductor. The external 
magnetic field H is directed along the c axis. We as- 
sume that H c i < H <C H C 2- In these conditions, the 
Abrikosov vortices are well separated and the amplitude 
of the superconducting order parameter is constant al- 
most everywhere, except from the core regions, whose 
size is of the order of the coherence length, £o ~ 1 nm. 
The average distance between vortices, a v , is estimated 
from the requirement that there is one flux quantum per 
vortex: a v = ^JWqJttH , where $o = hc/2e. In typical 
experimental situations, £o <C a v <C Xl, where Xl is the 
London penetration depth, which allowed us to replace 
the internal induction B by the applied field H . 

The quantity measured in STM experiments is the lo- 
cal differential tunneling conductance, which is propor- 
tional to the local DoS, N(r,u) = -(1/tt) Im Gf^r, r; lu), 
where G R is the retarded Gor'kov-Nambu matrix Green 
function. In the presence of a single scalar impu- 
rity, one can express G R in terms of the Green func- 
tion, G R (ri, Vi\ lu), of a clean superconductor in the 
mixed state: G R = G R + G R TG R , where T(w) = 
UT3 [1 — ugo(oj)T^\ 1 is the T-matrix, with go (u>) = 
G R {0,0;w). Then, 



N(t,u) = N (r,u}) + N imp (r,w), 



(1) 



where No is the local DoS for a d-wave superconductor 
in the mixed state in the absence of impurity, and 

N imp = -\ Im [G R (r, 0; w)T(w)Gf (0, r; lu)] n (2) 

is the impurity-induced contribution. The unitary limit 
corresponds to strong scattering: u — > 00 and T(uj) — > 
— g^ 1 ^). In this limit, iV(0,cj) = because the quasi- 
particles are prevented from occupying the impurity site 
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by a strong repulsive potential. Thus, in order to study 
the effect of magnetic field on the bound state near a 
unitary impurity, one should calculate the local DoS 
at one of its nearest neighbors, e.g. at r = a, where 
N(r, uj) reaches its maximum That the STM pictures 
show the maximum of the tunneling conductance directly 
above the impurity can be attributed to the blocking ef- 
fect of Bi-0 layers 0] and does not invalidate the stan- 
dard theoretical model of the impurity bound states. It 
should be mentioned that there are alternative explana- 
tions of this effect. For example, it was suggested that 
superconductivity is completely destroyed in the vicin- 
ity of a Zn impurity (the "Swiss Cheese" model, see Ref. 

creating an effectively normal region with high DoS. 
Another possibility is the Kondo screening of the local 
magnetic moments induced around the Zn impurity, see 
Ref. [l6|. It is still an open question, whether the peak 
in the STM conductance right above the impurity site is 
related to the details of the tunneling measurements or is 
an intrinsic property of the Cu-0 layers. In this article, 
we adopt the former point of view. 

We shall see that the energy dependence of N imp (a, uj) 
is determined mostly by that of the T-matrix. A closed 
analytical expression for go(io) = Gq '(0,0; w) can be ob- 
tained only in the case of a uniform supercurrent, when 
a gauge transformation can be applied to make the order 
parameter real and independent of r, so that the Fourier 
transform of G§(r, 0; to) is given by 



where z = toj Aq and 



(uj + - v F q s )r + £kT3 + A k n 



VFq s 



(3) 



Here w+ = lu + iO, r< are the Pauli matrices, £k is the 
normal state spectrum, Vp = Vk£k is the Fermi velocity, 
q s is the superfluid momentum, and Ak = 2A (cos k x d— 
cos k v d) is the mean-field order parameter, corresponding 
to d x 2_ y 2 symmetry (d is the lattice constant). We do 
not calculate the order parameter self-consistcntly and 
assume Ao to be constant. The numerical investigation 
of the self-consistency effects shows some suppression of 
the order parameter near the impurity site |l7j |. which 
leads only to a renormalization of the effective impurity 
strength towards the unitary limit [l^. It follows from 
Eq. JSJ that the quasiparticle energy is Doppler-shifted in 
the presence of supercurrent: E^ = ±\/£k + A^. + Vi?q s . 

The relevant energies of the problem are small com- 
pared to the gap magnitude: according to Refs. flflj |. 
the energy scale associated with the Doppler shift in the 
mixed state is given by Eh = vp/a v <C Ao. The Zeeman 
splitting is neglected because fJ-sH <C Eh at all experi- 
mentally relevant fields along the c axis (the field at which 
the Zeeman splitting becomes comparable to Eh is of the 
order of Hz jEL c i , where H p is the paramagnetic limiting 
field). Assuming an electron-hole symmetric band, the 
momentum integrals can be easily calculated, giving the 
following result at real frequencies: go(uj) = ttNfF(z)tq, 



m = E 



1 i 
— {z - z n )\n\z - z n \ - ~\z - z n \ 
Zn 4 



(4) 



Here Nf is the DoS in the normal state at the Fermi 
level, n labels the gap nodes, z n = v n q s /Ao are the 
dimensionless Doppler shifts, and v„ are the Fermi ve- 
locities at the nodes (we label the nodes in such a way 
that Z\ = — Z3, and z 2 = —24). Expression Q is valid at 
\z\ < 1, \z n \ < 1. Finally, T(z) = uct z [c - F{z)t z ]- 1 , 
where c = 1/(ttuNf) = cot<5o > for a repulsive impu- 
rity. 

The spectrum of the bound states is determined by the 
poles of the T-matrix, i.e. by the equation F(z) = ±c. 
In zero field, F(z) -> F (z) = (2/n)z1n\z\ - i\z\. At 
c 1, the equation Fq(z) — ±c has the solution 
zq = =F7rc/(2| lnc|) — i7r 2 c/(4 In 2 c) with logarithmic ac- 
curacy, describing a narrow impurity-induced resonance 
In the unitary limit c = 0, the resonance is decou- 
pled from the continuum of bulk excitations and becomes 
a true bound state of zero width, which manifests itself 
as a sharp peak in the tunneling DoS. In the presence 
of a non-zero supercurrent, the DoS of the bulk excita- 
tions does not vanish at lu — 0, which leads to a stronger 
hybridization of the bound state and the continuum of 
propagating states. As a result, the bound state is re- 
placed by a sharp resonance whose width is proportional 

to v F q s /\ln(v F q s /Ao)\ < A 0. 

In the mixed state, the supercurrent and the order pa- 
rameter amplitude are non-uniform, and Eq. @ is no 
longer valid. Further progress can be achieved by using 
the so-called Doppler-shift approximation 0, 0] . The 
basic assumption of this approach that the low-energy 
properties of a superconductor with the gap nodes are 
determined by the extended quasiparticle states and the 
contribution from the vortex cores can be neglected. The 
effect of the magnetic field can be described quasiclassi- 
cally by introducing a non-uniform Doppler shift of the 
excitation energy: E^ — > -Ek(r) = £4. + Vi?q s (r). In 
particular, the Green function at coinciding arguments 
becomes Gq{y,v;lo) = itNfF(z)t , where F is defined 
by Eq. Q, in which the Doppler shifts z n now depend 
on the coordinate: z n — > z n (r) = v„q s (r)/Ao. It is ex- 
pected that the Doppler-shift approach works well at low 
fields H -C H C 2 (for a discussion of its applicability and 
limitations, see, e.g., Ref. [20|). 

The Doppler-shift approximation should be used with 
caution in the calculation of the local DoS because 
the impurity-induced contribution rapidly oscillates as 
a function of r with the period of the order of A;^ 1 . We 
assume that the impurity is unitary, so that 



Im 



[G£(r,0;w)G*(0,r;w)] u }. (5) 
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Due to the smallness of the energy scale Eh, only 
the limit (uj,vf<1s) < Eh <C Ao is relevant. In this 
limit, one can neglect Nq(t,u) compared to N imp (r,uj), 
which is peaked at small u) because of the singularity of 
F~ l (u)/Ao). In addition, the product of two Green func- 
tions in the second term on the right-hand side of Eq. 
iJSJ depends weakly on lo and q s and can therefore be 
replaced by its value at lo = q s = 0: 



[<(r,0;0)G«(0,r;0)] n = (vrJV F ) 2 7 (r). 



(0) 



Here 7 = |Ji| 2 + \ I2\ 2 is a real number, which depends on 
the order parameter and the shape of the Fermi surface, 



1 r d 2 k 

ttN f J (2tt) 2 
1 f d 2 k 



6, 



ttN f J (2tt) 



where the k-integration goes over the first Brillouin zone 
—ir/d < k x , k y < n/d. Thus, in the unitary limit and at 
w < A , 



A^ mp (a,u) 
N F 



7 (a) Imr 1 ( — 



(7) 



For a dispersion £k = — 2t (cos k x d + cos k y d) — 
At\ cos k x d cos kyd — /x, and r = a = xd, we have 7(a) = 
0.63 (we used the parameters representative of BSCCO 
near the optimum doping: t\jt = —0.3, (j,/t = —1, 
A /< = 0.04, A = 10 meV). 

When calculating the energy dependence of the impu- 
rity resonance contribution to the local DoS, Ni mp (uj) = 
Ni mp (a.,uj), one can use the values of the Doppler shifts 
at the impurity site, z n = z n (0), since the variation of 
the supercurrent at the atomic length scale is negligibly 
small. As a final step, one should average over all impu- 
rity positions (or, equivalcntly, over the vortex positions 
relative to the impurity site): 



N lmp (uj) 
N F 



(8) 



Here (ImF" 1 )? = / dz 1 dz 2 V(zi, z 2 ) ImF" 1 , where F 
is defined in Eq. Ijljl. and V is the distribution function 
of the Doppler shifts: 



v 2 q s (r) 



xS ( z 2 



An 



(9) 



where A is the system area. 

According to Ref. Zn impurities in BSCCO act 

as pinning centers for the vortices, so that a sizeable 
fraction (~ 50%) of the vortices reside at the impurity 
sites, creating a disordered vortex solid. Because of a 
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FIG. 1: The suppression of the impurity contribution to the 
local DoS at the nearest-neighbor site r = a, at different 
values of the magnetic field: zh =0.1 (solid line), zh = 0.2 
(dashed line), and zh = 0.4 (dotted line). 



strong variation of the order parameter in the vortex core, 
the Doppler-shift approach is clearly inapplicable for the 
bound states near such impurities. On the other hand, 
in the experimental conditions of Ref. ^(Ji the majority 
of Zn atoms are found in the space between the vortex 
cores and exhibit zero-bias peaks in the tunneling DoS. 
The field dependence of these peaks can be analyzed us- 
ing the Doppler-shift approximation. The overlap of the 
impurity resonances and the core states is neglected. 

The probability distribution, Eq. {SJ, of the Doppler 
shifts at the inequivalent gap nodes depends on the vor- 
tex configuration. In the case of a disordered vortex solid, 
which is relevant to the experiment, the probability den- 
sity is peaked around the zero Doppler shift. In Ref. [20| . 
different models for the supercurrent distribution were 
discussed giving rise to different V{z\,z<i). One possible 
choice, which agrees with the results of numerical simu- 
lations and is believed to capture the physics reasonably 
well, is given by 



1 



T( Zl ,z 2 ) = - 



7T (z\ + Z% + Z|) 2 ' 



(10) 



where zh — Eh/^q — \J H/ H c2 . We use this distribu- 
tion in Eq. ||SJ to calculate the average DoS. 

At zero field, Zh — 0, and V(zi,z 2 ) — 8(z\)5(z 2 ). In 
this case, the local DoS has a divergent peak near zero 
bias: Ni mp {Lo)/Np = jIkiFq (z) [it is not a <5-function 
peak because F ~ 1 (z) does not have a simple pole struc- 
ture]. At finite H, the height of the peak at u> = 
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becomes finite and can be calculated exactly: 



N imp (0) 
N F 



2l 



1 



pi 



22 



2V21n(V2 + l) 7 /Sc2 

£ ^VT 



For the band dispersion described above, N imp (0)/Np ~ 
1.57 t/H&Jh. This should be contrasted to the uniform 
supercurrent model considered in Ref . , according to 
which N imp (0)/N F oc fT -1 . 

It can also be shown that the result Ni mp (0)/Np oc 



y/H^/H is actually independent of the choice of 
"P(£i, 22)- Indeed, "P(zi, 22) is even in both zi and Z2 and 
also symmetric with respect to the interchange z\ ^ z 2 . 
If the only energy scale of the Doppler shift distribution 
is Eh, then V can always be written in the form 



V(z 1 ,z 2 ) = — / — , — 
z H \z H z H 



(12) 



where f(x,y) is a dimensionless function such that 

ffav) = f(y> x )> f{ ±x ,±y) = f(x,y), and 
J_ oc dx dy f(x,y) = 1 in order to satisfy the normaliza- 
tion condition for V . Then, 



N imp (0) 



27 

z 2 



dz\ dz2 p f z i 



Z2 



f 



£2 

\zh' z H 
Ao 



A- 



E 



H 



(13) 



where A = 27 J da;dy/(a;,y)/(|a;| + 

At arbitrary bias, the impurity contribution to the 
local DoS given by Eq. © can be calculated numeri- 
cally. In Fig. ^ we pl°t N imp (uj)/NF for several val- 
ues of zh- The most prominent feature of these graphs 
is that the peaks are indeed quickly suppressed as H 
grows, in accordance with the result 1)11(1 . To see that 
this suppression is accompanied by broadening and to 
facilitate comparison of the results at different zh, we 
plot in Fig. |21 the local DoS normalized to its maxi- 
mum value, Ni mp (u)) / Ni mp (0) . The broadening of the 
peaks can be attributed to the hybridization of the im- 
purity bound states with the field-induced nodal quasi- 
particles, whose density of states at zero energy increases 
with field as yj H/ H C 2 [l2( . One should keep in mind that 
since the expressions (@J and JSJ are applicable only at 
z, zh <§C 1, they give a quantitatively correct picture only 
in the vicinity of the zero-bias peaks at fields far from 

EI C 2- 

In conclusion, we have studied the influence of 
magnetic field on the impurity bound states in d-wave 
superconductors. The magnetic field effects are treated 
semi-classically, using the Doppler-shift approach. The 
main result is that the bound states are destroyed by 
magnetic field, which manifests itself as the suppression 
and broadening of the zero-bias peaks in the tunneling 
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FIG. 2: The broadening of the DoS peak at r = a caused 
by magnetic field. The normalized average DoS is calculated 
at different magnetic fields: zh = 0.01 (bottom curve), zh = 
0.02, z H = 0.05, z H = 0.1, z H = 0.2, z H = 0.4 (top curve). 
Since the peaks are symmetric with respect to u — > — ui, only 
the positive bias is shown. 



DoS. The height of the peaks scales as H^ 1 ^ 2 for any 
non-uniform supercurrent distribution, depending on a 
single energy scale Eh- This result is sensitive to the 
orientation of the field: if H is parallel to the ab planes, 
then the uniform supercurrent model of Ref. |l3| should 
be more appropriate, which predicts that the height of 
the peaks is proportional to H . As mentioned above, 
the impurity resonance model is not the only possible 
explanation of the zero-bias anomalies in the tunneling 
conductance of HTSC's. It would be interesting to 
look at the magnetic-field response in the other models, 
which might help to resolve the controversy about the 
origin of the zero-bias peaks. 
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